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ABSTRACT 

We determine the slope of the power spectrum of the matter density perturba- 
tions from the Tenerife observations of the cosmic background radiation temperature 
anisotropies. We compute the projected radiation anisotropy power spectrum measured 
by this experiment and study its dependence with respect to the slope of the temperature 
anisotropy power spectrum on the sky. We show that Tenerife alone implies the upper 
bound on the spectral index of m < 3. Stronger conclusions can not be reached due to 
the small data set. The method proposed can be applied to any small scale experiment. 
Sampling the same region of the sky with different window functions could probe the 
slope of the radiation anisotropy power spectrum at different scales and confirm the 
presence of Doppler peaks. 

Subject headings: Cosmic Microwave Background. Cosmology: theory. Methods: nu- 
merical 



1 



1. Introduction. 

In the last few years several groups have reported 
measurements of temperature anisotropies on the 
Cosmic Microwave Background (CBR) at different 
angular scales (see White, Scott & Silk 1994; Bond 
1997, for a review). These measurements are be- 
coming sensitive enough to disproof models of large 
scale structure and to measure cosmological parame- 
ters such as the amplitude and slope of matter den- 
sity perturbations (Bennett et al. 1996), the height 
and location of the first Doppler peak in the temper- 
ature anisotropy power spectrum (Scott et al. 1995), 
and to determine the possible existence of a cosmo- 
logical constant (Lineweaver et al. 1996). Assum- 
ing temperature fluctuations of the CBR to be Gaus- 
sian distributed, the power spectrum of the tempera- 
ture anisotropies contains all the physical information 
about the temperature field. The power spectrum has 
been obtained for the COBE/DMR full sky temper- 
ature anisotropy map, where window function facil- 
itates the decomposition of the temperature fleld in 
its spherical coefficients (Bunn, Hoffman & Silk 1996). 
The variance of data from small scale observations has 
been used to constrain models of Large Scale Struc- 
ture (Ganga et al. 1996; Gorski et al. 1996; White et 
al. 1996). 

To the best of our knowledge, no attempt has been 
made to estimate the slope of the radiation anisotropy 
power spectra for experiments with partial sky cover- 
age but higher resolution than COBE/DMR. In this 
article we shall perform such analysis aiming to probe 
the overall shape of the power spectrum at intermedi- 
ate and small angular scales. We shall center our dis- 
cussion on the analysis of one-dimensional data sets 
and, in particular, on the scan measured by the Tener- 
ife group (Hancock ct al. 1994). The situation is com- 
mon in cosmology where information about the mat- 
ter power spectrum is inferred from one-dimensional 
pencil beam galaxy surveys (Broadhurst et al. 1990) 
or from the angular correlation function of galaxies 
(Maddox et al. 1990). The outline of the paper is 
as follows: in Section 2, we show that for small scale 
experiments the power contained in the data is re- 
lated to the temperature anisotropy power spectrum 
on the celestial sphere. In Section 3 we will assume 
that the radiation anisotropy power spectrum is given 
by a power law with slope m, roughly constant on the 
scales probed by the Tenerife experiment. We shall 
show how this slope can be measured from the data. 



In Section 4 we apply the technique devised in the pre- 
ceeding section to the temperature anisotropies ob- 
served at Tenerife. We perform several Monte Carlo 
simulations to determine the confidence level of our 
estimates. Finally, in Section 5 we discuss our results 
and present our main conclusions. 

2. Power Spectrum of one-dimensional tem- 
perature anisotropy fields. 

If temperature anisotropies on the sky arc Gaus- 
sian distributed, they are fully described by the an- 
gular correlation function 

C{6) = j^J2^2l + l)CiWiPi{cos6), (1) 

where C/ = |a/mp is the radiation anisotropy power 
spectrum and Wi is the window function for a par- 
ticular experiment. In the limit / ^ 2 the correla- 
tion function can be rewritten as (Bond & Efstathiou 
1987): 

C{0) = ^ I ldlC{l)W{l)Jo{ie). (2) 
2-K J 

Jo denotes the cylindrical Bessel function, and C{l),W{l) 
represent continuos functions that verify C{1) = Ci,W{l) = 
Wi at positive integers. This limit is equivalent to 
consider the sky as a flat surface and within the same 
approximation the temperature anisotropy can be ex- 
pressed as 

where r is a direction on the sky. The relation between 

the " continuos radiation anisotropy power spectrum" 
defined by Prad{k) = \ST/Taik)\'^ and C{1) is 

k^Pradik) = PC{l)W{l). (4) 

In this last expression, angular and linear scales arc 
related by r = RhOJ = kRn, being Rh the horizon 
radius that, for convenience, we shall take equal to 
unity. 

Most observations on small angular scales sample a 
small fraction of the sky. In general, the experiments 
measure the temperature along a single angular direc- 
tion. In the flat sky approximation the temperature 
anisotropy along, for example, the x axis is 
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All the statistical information contained in the data is 
determined by the one-dimensional power spectrum 
Pioiq) = |(5T/To(g)|2, being ST/To{q) the ID Fourier 
transform of the data. Since the correlation function 
is common to both the 2D and ID temperature field 
because of statistical isotropy, the projected power 
along one direction and the sky power spectrum are 
related by the integral equation (Peacock 1992) 

PM^) = ll dk2Prad{\Jq'+kl). (6) 

In our approximation, equation (6) can also be ex- 
pressed as 



(7) 



l>q 



The last identity follows from equation (4). An im- 
mediate consequence of this equation is that P\D{q) 
must be a decreasing function of wave number. 

In principle, inversion of equation (7) would give 
the radiation anisotropy power spectrum on the ce- 
lestial sphere C{1). This type of inverse problems are 
common in astronomy and are usually ill-conditioned 
(Lucy 1994). The number of data points is usually in- 
sufficient to recover the power. Furthermore, observa- 
tional errors in the data are amplified by any inversion 
procedure and the final answer is totally unreliable. 
Instead of trying to recover the full power spectrum, 
we will assume it to be described by a power law at 
the scales where the window function is non-negligible 
and we shall discuss in the next section how the slope 
can be recovered from the data. 

3. The slope of radiation anisotropy power 
spectrum at Tenerife scales. 

As mentioned in the introduction, we will center 
our study on the Tenerife experiment. For most cos- 
mological models with fio = 1, temperature anisotropies 
at Tenerife scales are dominated by fiuctuations in 
the gravitational potential (Sachs & Wolfe 1967). For 
matter spectral indexes n < 3 the relation is given by 
(Peebles 1982) 



Ci = (aa) 



l)/2]r[(9-n)/2] 



r[Z + (5-n)/2]r[(3 + n)/2]' 



(8) 



The quadrupole 02 will be taken here as a normaliza- 
tion constant. The radiation anisotropy power spec- 
trum in the region 2 << Z < 50 is essentially a power 



law: P{k) ~ fc™ and its slope m is related to the 
matter spectral index n by: m « n — 3. 

One can estimate quantitatively the dependence 
of the radiation anisotropy power spectrum PiD{q) 
with m on scales 20 < Z < 40. In that region, the 
power spectrum is Prad{k) « fc™ exp (— fc^cr^) because 
the window function is dominated by the exponential 
cut-off introduced by the finite beam size (a) of the 
antenna. Within this approximation the projected ID 
power spectrum defined by (eq. [6]) is 



1 

PiD{q) = -e-'^^l dy{y'+q'y 

Jo 



^e-y . (9) 



This integral can be calculated in closed form in terms 
of Whittaker functions (Gradshteyn & Ryzhik 1980): 



P,n{q) = ^^e(-i-^«^)gi('"-i)Ty,,_,(gV2), 

(10) 

with a = {1 + m)/4. In the limit qa > 1, the depen- 
dence of the power on q reduces to: 



Pic(9) ocg"exp(-a2g2). 



(11) 



This is our main result of this section: at high multi- 
poles the projected power spectrum Pioiq) depends 
on the spectral index in the same way as the radiation 
anisotropy power spectrum Prad{k). The ID power 
is boosted by the projection of short wavelengths in 
two dimensions, like in pencil beam surveys (Kaiser 
& Peacock 1991), but the overall shape remains un- 
changed. This result should not come as a surprise 
since it is a consequence of the statistical isotropy of 
the correlation function. Nevertheless, it has an in- 
teresting application to observations. On the region 
dominated by the exponential tail of the window func- 
tion the amount of information contained in a 2D map 
or in a scan is similar. Therefore, experiments like 
Tenerife or Bartol-IAC (Piccirillo et al. 1996) that 
obtain data in the form of a single scan can constrain 
the radiation anisotropy power spectrum at the same 
confidence level than maps do. The main limitation 
does not come from the geometry of the experiment, 
but by the number of independently measured data 
points. 

Since we will analyze the data measured by the 
Tenerife group, we will restrict our analysis to this 
experiment. The Tenerife group is continuously mea- 
suring a strip of 360 degrees in R.A., at declination 
40°. Their observational strategy is a double beam 



3 



switch with an 8° chop in right ascension with 5.5° 
beam width. The window function of the experiment 
peaks at Z = 20 and is rather insensitive to multipoles 
outside the range 10 — 50. Their final data is displayed 
as a scan of fixed declination and right ascension from 
12 to 18h. Outside this range, the signal is dominated 
by the galactic contribution and other systematic er- 
rors. For the Tenerife experiment we shall use infor- 
mation contained in the range 20 < Z < 40. For small 
multipoles the curvature of the sky renders our treat- 
ment of the data invalid. For multipoles larger than 
I = 40 the window function erases the cosmic sig- 
nal and the data contains only uncorrelated noise. In 
the region of interest, \og{PiD{q)) is almost a straight 
line. In Figure 1 we plot the projected power spec- 
trum for different slopes: m = —2.5, —2, —1.5, 1. The 
corresponding spectral index of matter power spec- 
trum is approximately n ~ m + 3. The ID power 
spectrum was computed using equations (7) and (8) 
for m < and equation (6) for m > 0. The slope 
miD = Alog(Pic(g))/Ag was obtained by means of 
a least square fit to the one dimensional power spec- 
trum in the region I = 20 — 40. More than 50% of the 
total power measured by the experiment is contained 
in this region of I space. In column 2 of Table 1 we 
quote the figures for different slopes. We verified that 
the numerical value and the analytical estimate using 
equation (11) were very similar. 

It is not surprising to find that miu depends 
weakly on m. At high multipoles, the power is domi- 
nated by the finite beam width and the signal is expo- 
nentially damped. Even though, we shall see in the 
next section that meaningful statistical information 
can be obtained from the data. 

4. Estimation of the radiation spectral index 
from Tenerife data. 

Our procedure to estimate the slope of the mat- 
ter power spectrum measured by the Tenerife experi- 
ment is very simple. We compute the ID power spec- 
trum from the data by means of a Fourier transform. 
The slope miD is obtained by a least square fit to 
\ogPir>{q) in the region 20 < Z < 40. We compare 
this value with the theoretically expected value that 
can be read from Table 1. 

To estimate the statistical uncertainty associated 
to the numbers given in Table 1, we perform Monte 
Carlo simulations by selecting five hundred Tenerife 
scans randomly placed on a single realization of the 



sky for any given n = m + 3. We computed the ID 
power in order to estimate mm- Each scan consisted 
of 90 points separated by 1° in R.A. This data was 
later binned to obtain an independent point per beam 
width. Our simulations were run considering zero 
noise and were repeated by adding noise levels of 20, 
50 and 100% of the cosmic signal to each scan. The 
slopes were Gaussian distributed around the mean, 
given in Table 1. The first column gives the slope of 
the radiation anisotropy power spectrum. In the sec- 
ond and third columns we write the theoretical ex- 
pected value of mm and the value for our particular 
realization of the sky. On columns 4 to 7 we give the 
measured slope for different noise levels. One should 
notice that the effect of white noise in our simula- 
tions is to increase the slope (to decrease it in ab- 
solute value). This is easily understood since white 
noise has a slope m = and with increasing ratio 
(Tn/cTs the scan has a larger contribution from larger 
slopes. 

The shortcomings of power spectrum analysis based 
on periodograms or other types of spectrum estima- 
tor are well known (Press et al. 1992). For example, 
according to equation (7) Pioiq) must be a decreas- 
ing function of q and in general it was not. Noise and 
data discretness and ahasing between different scales 
induce oscillations in the one dimensional power that 
results in a very high (of the order of 0.8) r.m.s. de- 
viation {a) on the values of columns 4-7 of Table 1. 
With this result, it would be impossible to distinguish 
between m = —2 and m = 1 models at la confidence 
level, even for a noise free scan. A partial solution to 
the problem can be obtained by averaging power spec- 
tra from several data sets (see also Press et al. 1992). 
Therefore, we repeated our simulations but averaging 
the power of N different scans before computing the 
slope. We took N' = 3, 5, 10 and the result was that 
the mean slope of the 500 scans did not change and 
the variance was reduced by a factor proportional to 
■s/TV. The values given in Table 1 were obtained after 
averaging three scans before computing the power. 
The r.m.s. deviation was a « 0.04 and was never 
larger than 0.05 in all cases. 

We finally apply our formalism to the data. The 
Tenerife group has pubHshed the data on three dif- 
ferent frequencies at declination 40° and soon will be 
releasing their measurements for two more declina- 
tions. The data at dec. 40° is displayed in Figure 
2. We plot the 10 GHz, 15 GHz and 33 GHz scans 
from top to bottom. In the same figure, the smooth 
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curves correspond to a Wiener filtered data set (Press 
et al. 1992). Our approach, based on calculating the 
power spectrum of the original data, allows for the 
use of filtering techniques to remove noise. Given a 
prior model and a noise level, the Wiener filter gives 
an optimal estimate of the background signal. In Fig- 
ure 2, the Wiener filtered data was constructed using 
a n = 1 prior model and a noise level of 10/iK, ap- 
proximately the noise level quoted by Hancock et al. 
(1994) for the 33 GHz scan. In Figure 3 we plot the ID 
power spectrum obtained from the Fourier transform 
of the raw data (upper solid curve) and the Wiener 
filtered data (lower solid curve) . The data was binned 
before computing the power, and the spectrum of the 
three frequencies was averaged to compare with the 
results of Table 1. For comparison, in Figure 3 we 
also plot, with dashed lines, the theoretical power for 
m = -2, 1 (n = 1,4). 

One should notice the high frequency tail present 
in the raw data. This high oscillations are due to the 
noise present in the original data. They do not have 
a cosmological origin since above I = 50 the window 
function goes to zero. However, this high frequency 
tail drives the slope to higher values. To see that 
this is indeed the case, we estimated the slope for the 
Wiener filtered data. The effect of the filtering was to 
fully remove the high frequency tail and to bring the 
slope close to to = —2. In Table 2 we give the slope 
of the ID power spectrum obtained from the original 
Tenerife data for the three different scans of 10, 15 
and 33 GHz and from the Wiener filtered data. We 
should remark that the raw data is compatible with 
high values of to indicating a non-negligible noise con- 
tribution to the temperature anisotropy. For compar- 
ison, assuming a prior model of n = 1 and a noise level 
of about 25%, the slope of the power spectrum rises 
significantly for the scan of 15 GHz and comes close 
to a value compatible with to = — 2 (n = 1). The 
average of the 3 scans is, in the raw data, compatible 
with TO = — 2, even though it centers on values around 
TO = 0, indicating that white noise is present in the 
data. In Table 2 the error bar attached to the slope 
averaged over three scans comes from the result of the 
previous Monte Carlo simulations. One can see that 
from the Tenerife data alone, a cosmological signal of 
TO > 3 can be ruled out at la level. On the other 
hand, for the Wiener filtered data, mio is very close 
to TO = —2. Furthermore, at Icr level, the Wiener fil- 
tered data is incompatible with to = 0. If the Wiener 
data is an optimal reconstruction of the signal mea- 



sured by the experiment, this result argues towards a 
cosmological origin of the temperature anisotropy in 
the scan. 

5. Discussion. 

In this article we have seen that on small angular 
experiments information can be obtained about the 
slope of temperature anisotropy power spectrum. It 
is clear, from the results of our previous section, that 
our estimates are less conclusive than COBE/DMR. 
However, this is not due to analyzing a ID scan in- 
stead of a full 2D map but it is due to the lack of a 
large number of independent data points in the Tener- 
ife data, about 20 per frequency channel. Our estima- 
tion of the radiation anisotropy power spectrum was 
done independent of other experiments and in a dif- 
ferent region of Z-space than COBE/DMR. Our main 
conclusion is that the slope is to < 3 at Icr confidence 
level. 

Experiments like Max and MAXIMA (Hanany et 
al. 1996) or the combined observations of Tenerife and 
Bartol-IAC, that sample the same region of the sky 
but with different resolution, could benefit from an 
analysis like the one done here. Improvement in the 
quality of the data and larger data sets would decrease 
the uncertainties and to improve our estimation of the 
spectral index. Moreover, the temperature anisotropy 
power spectrum changes slope very quickly close to 
the Doppler peaks. For a CDM model with fij/i^ = 
0.02 the slope rises from m = — 2 on large angular 
scales up to m = — 1 before decreasing with to = — 3 
after the first Doppler peak. Experiments carried out 
on the same region of the sky, with high signal to 
noise ratio, at different frequencies, can be used to 
probe the existence of a Doppler peak. 
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scan 


data 


Wiener filtered data 


10 GHz 


- 0.070 


- 0.108 


15 GHz 


- 0.101 


- 0.185 


33 GHz 


- 0.148 


- 0.170 


10+15+33 


-0.11 + 0.05 


-0.15 + 0.05 



Table 2: Slope of the ID power spectrum for the three 
Tenerife scans. Column 2 corresponds to the original 
data and Column 3 corresponds to the Wiener filtered 
data where the prior model is n = 1 and the assumed 
noise level was 10/iK. 



m 


miDitheo) 


miDisim) 


^noise — 


^noise — 0.2(Ts 




^noise — 


-3.0 


-0.178 


-0.183 


-0.177 


-0.166 


-0.133 


-0.081 


-2.5 


-0.165 


-0.169 


-0.171 


-0.162 


-0.133 


-0.086 


-2.0 


-0.152 


-0.156 


-0.160 


-0.156 


-0.130 


-0.089 


-1.5 


-0.139 


-0.143 


-0.152 


-0.148 


-0.128 


-0.090 


-1.0 


-0.127 


-0.131 


-0.141 


-0.138 


-0.122 


-0.089 


0.0 


-0.104 


-0.104 


-0.119 


-0.117 


-0.106 


-0.082 


1.0 


-0.083 


-0.086 


-0.094 


-0.094 


-0.087 


-0.070 


3.0 


-0.047 


-0.050 


-0.052 


-0.052 


-0.049 


-0.041 



Table 1: In the second column we give the theoretical value of the slope of the ID power spectrum Pin{q) for the 
different radiation anisotropy power spectra. In the third column we give the slope of the single realization of the 
sky used to generate the scans. Columns 4 to 7 give the mean slope of five hundred simulated Tenerife scans. The 
values were obtained after averaging the power spectrum of three different scans before computing the slope. The 
r.m.s. deviation was very similar to all simulations {a « 0.04 — 0.05). 




Fig. 1. — One-dimensional power spectra (eq. [6]) as it would be measured from a Tenerife scan. Four different 
slopes have been considered. The power spectra displayed from left to right corresponds to: m = —2.5, —2, —1.5, 1. 
The most negative m gives rise to the largest negative slope. 
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16 14 12 10 

Right Ascension 



Fig. 2. — Temperature anisotropies of Tenerife scans at 10, 15 and 33 GHz. The smooth curve corresponds to the 
Wiener filtered data with a n = 1 prior model and noise level of 10/iK. 
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Fig. 3. — The solid lines correspond to the average power spectrum obtained from the three scans of original data 
(upper solid line) and of Wiener filtered data (lower solid line) . Notice the long tail of high frequency noise present 
in the original data. For comparison, the theoretical one-dimensional power for n = 1,4 (left, right dashed curve) 
are also plotted. 
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